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General local spin S ground states, described by a Valence Bond Solid (VBS) on a two dimensional 
lattice are studied. The norm of these ground states is mapped to a classical 0(3) model on the same 
lattice. Using this quantum-to-classical mapping we obtain the partial density matrix pa associated 
with a subsystem A of the original ground state. We show that the entanglement spectrum of p\ 
in a translation invariant lattice is related with the spectrum of a quantum XXX Heisenberg model 
and all its conserved charges on the boundary of the region A. 
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I. INTRODUCTION 

Quantum entanglement, the spooky action at a dis- 
tance that has been signaled as the characteristic of 
quantum mechanics^, has received renewed attention re- 
cently, specially with the growth of quantum information 
science^, and as a new tool to study properties of many- 
body systems^} It has been found that entanglement is 
sensitive to topologically ordered states^, quantum phase 
transitions^, and even magnetic properties of solids^. 

A complete description of the entanglement proper- 
ties of a bipartite pure state I*) composed of subsys- 
tems A and B is given by the entanglement spectrum 
(ES)P, i.e. the eigenvalues of the reduced density ma- 
trix (RDM) of subsystem A (or B). The RDM for the 
subsystem A, pa, is obtained by tracing out the de- 
grees of freedom belonging to B from the density matrix 
p = \^>){^>\ which characterizes the state |\E'). In general 
the density matrix of the subsystem A may be written as 
Pa = exp(— (3H e g), where H c s is an effective (also called 
entanglement) Hamiltonian. It has been shown that the 
entanglement Hamiltonian describes excitations living at 
the edge of partitions of the ground state of fractional 
quantum Hall states^, one dimensionaP^ and topolog- 
ical systemiP^EI In two dimensions, few analytical re- 
sults on the ent angle ment spectrum in generic spin sys- 
tems are known^ED Numerical studies of the ES have 
been performed in the 2D Affleck, Lieb, Kennedy and 
Tasaki (AKLT) modeP, which possess a known valence 
bond solid (VBS) ground state. 

Lou et aP3 and Cirac et aZP^l showed that the ES of 
a partition in the ground state of the AKLT model is 
related with the conformal XXX Heisenberg model on 
the boundary of the partition by using Montecarlo and 
projected entangled pair states (PEPS) in finite size sys- 
tems. In this paper we show that the ES of a partition 
of a whole class of ground states defined in translational 
invariant lattices, can be approximated by the thermal 
spectrum of a scries of local Hamiltonians which are the 
conserved charges associated with the XXX Hamiltonian 
defined on the boundary of the partition. 

In this paper we introduce the spin S model in sec- 



tion [TTJ defined on a two dimensional lattice wrapped 
on a tor us an d construct its explicit VBS ground state 
followin g ^ 20 * 23 !. Then, we derive an expression for the 
RDM (also called partial density matrix) pa in section 
III| This operator is expressed in terms of classical vari- 



ables in section IV In this representation, the operator 



can be expanded in different graph contributions of the 
classical 0(3) model as presented in section [V] From this 
expression we identify the Heisenberg Hamiltonian for 
spin 1/2 particles in the boundary as the leading term 
in a sequence of boundary Hamiltonians. In the last sec- 
tion, we summarize the results and discuss further possi- 
ble generalizations. 



II. 



SPIN S VBS GROUND STATE ON A TWO 
DIMENSIONAL TORUS 



As discussed orPHHHI] ^ j g p 0S sible to construct a va- 
lence bond solid (VBS) ground state in a planar graph 
Q (without loops) in the following way: Given a pla- 
nar graph Q, consisting of a set of vertices (sites) V and 
edges E, with Zi edges arriving to vertex i (in graph the- 
oretical language, Zj is called coordination number), we 
place a local spin Si on the vertex with the condition 
Si = Zi/2. The local spin state is constructed from the 
symmetric subspace of Zi spins 1/2 (doing this we obtain 
a higher spin representation of dimension 2Si + 1 from 
2Si fundamental representations of SU(2)). Finally we 
antisymmetrize between nearest neighbors. Represent- 
ing the spin 1/2 constituents of the spin Si at site i as 
black dots, using a circle to indicate symmetrization and 
a bond between antisymmetric neighbors, we obtain a 
planar graph Q' isomorphic to Q, see Fig[l] 

The AKLT Hamiltonian for which the VBS state con- 
structed is a ground state is a sum over interactions on 
all edges E of Q, H = J2{k,i)eE H ki{Sk + Si), where the 
Hamiltonian density H k i is 
s k +s l 

H kl (S k +Si)= A ki^k + S l ), (1) 

J=S k +Si + l-M lk 

the coefficients A'^ > are arbitrary and can depend 
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FIG. 1. a) Original planar graph Q, vertices represented with 
black dots, b) VBS state on Q' , circles (vertices of Q') rep- 
resent symmetrization of constituents spin 1/2 (black dots) 
particles, while bonds represent anti-symmetrization of neigh- 
boring spins. 



factor. \A a ) and \B a ) are states defined in the subsys- 
tems A and B respectively, with the total system being 
the union of A and B. This states {|^4 a )> |-B«)} are a 
complete, linear independent but not orthonormal ba- 
sis in principle. The density matrix for this pure state 
is the projector onto the ground state p = Af\^f)(^\, 
with M^ 1 = (^1^). Tracing out the sites belong- 
ing to the subsystem B, we obtain the partial den- 
sity matrix, which describe the system A, pa = Tr^p. 
Using pi) the partial density matrix becomes pa = 
NYl a p\Bp\Ba)\Aa) ® (Ap\. We can write the opera- 
tor pa as a matrix using the basis \i) — J2 a Ui a \A a ) and 
its dual (j\ = ^ 7 Vj~ f {A 1 \, where U and V are matrices 
of change of basis and 



Sij. We have (j\pA\i) = 
0>a)h = ATY Ja ^{ B p\ B c t )V jl {A 1 \A a )U ill {Ap\A li ). 
However from the orthonormality condition = 5ij 
follows the relation ^ 7 Vj 1 {A 1 \A a ) = (U^ 1 )^ which in- 
serted back in the expression for (pa) ji simplifies it to 
(j\PA\i) =^T, a ^U i ^{A^\Ap))*{Bp\B a ){U- l ) aj . As 
U is a matrix of change of basis, (pA)ji is just a similar- 
ity transformation of the operator (pA)fj,a given by 



(pa)^ = 



((A^Aj))* (B 7 \B a ) 



(4) 



on the edge (k,l), while the positive number M^i is the 
number of bonds (edges) connecting the sites k and I. 
The operator TTjy (Sk + Si) is a projector of the total spin 
Jki = Si + Sk of the edge (k, I) on the subspace of spin 
value J, its explicit form is 



hi) 



Sk+Si 

n 

j=|S fc -S 1 |,i#J 



(Jfcl) 2 -j(j + l) 
jki{Jki+i)-ju+iy 



(2) 



This matrix is not hermitian in the usual sense 0^ = 0, 
but it is isospectral with (pa) ■ Defining the overlap 
(Aa\Ayy), the partial density matrix 



matrices M, 



is constructed gluing together two overlap matrices, one 
for each subsystem, along the boundary of the partition, 
leaving one index free in each overlap matrix, obtaining 
a torus with a cut along the partition (see Fig. [2|. 

Using the Schwinger boson representation f or sp in op- 
erators, the VBS state on Q can be written a d 25 * 26 ! 



The VBS state is the unique ground state of B^. 
While this construction is totally general, in the rest of 
this discussion we will focus on graphs without bound- 
aries, which can be embedded on a two dimensional torus, 
with Mij = 1 for all edges. 



III. PARTIAL DENSITY MATRIX AND 
SCHWINGER BOSON REPRESENTATION OF 
VBS GROUND STATE 

Using the Schmidt decomposition, any ground state 
\^f) of a system can be written as 



I* 



|*} = ^|A Q )®|B a ), 



(3) 



vbs/ — 



= IT (a$-a$)|0>, 



(5) 



where Eg is the set of all edges (bonds) of Q and |vac) 
is the state annihilated by all the and hi operators, 
i.e. cij|0) = bi\0) = 0, V i. For a generic partition of the 
system into two subsystems A and B (we assume both of 
them to be simply connected) with boundaries dA and 
dB, we have a collection of vertices Va, Vb such that 
Va U Vb — Vg and a collection of edges (bonds) which 
endpoints live either both in A (B) or one in A and the 
other in B. For bonds which both endpoints live in A 
we will say £ Ea (similarly for B), while for shared 
bonds with endpoints i and j we use i £ dA, j £ dB. The 
set of shared bonds we will call it d (and is the same for 
A and B). Finally the cardinality of a set M is denoted 
by \M\. Using this definitions, we can write the state 
\^ vbs ) in the form Q 



where the states \A a ) and \B a ) are related to the usual 
states appearing in the Schmidt decomposition by a scale 



vbs / 



II (4b}-a}bl) JI (ajat + fot 6 t )|0>;(6) 
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3 



A 



V 



B 



FIG. 2. (color online) Graphical representation of the par- 
tial density matrix pA, for a particular partition. The overlap 
matrix Afjjl' = (A a \Ap) corresponds to a section of the torus 
with two different, inner and outer, boundaries. To compute 
(Pa)c/3, we glue the inner boundaries of the overlap matri- 
ces M' A ' and M' s ', obtaining a torus with a cut along the 
boundary of the partition. 



where we have applied a local basis transformation on 



the sites (vertices) in B, a] —> —b\ and fe - — ► at just for 
latter convenience. In the shared bonds, we can assign 
to an endpoint j of a bond, it's partner in the other end 
of the bond to be j. Doing this we can expand ^ in the 
forirPl 



\* vbs ) = £ n (<4r(bt) i - a <(<4r(bt) 1 - 

{a} iedA 

x n (at6j- a t 6 t)|0) 

(i,j)eE A UE B 
= £l^{a}>®|S {a} ), 



(7) 



{a} 



here {a} = {ai, a%, ■■oi\dA\\ with a, = 0,1, labels the 
different ground states of the subsystems, which span a 
Hilbert space of dimension 2' 9j4 l. The Hamiltonian in A 
is H A = J2jk,i)eE A H ki (and similarly for B), with H kl 
defined in ([I]) . From Q , we can read off the form of the 
states |^4{ Q }) 

\A {a} )= n n ( 6 i) i " a< («i) Qi io), (8) 

{i,j)eE A iedA 

using Q and we can compute the density matrix pa 
in terms of the overlap matrices Mr^l = (A{ a y\A{py). 
The computation of this overlap matrix can be mapped to 
the computation of partition and correlation functions in 
an 0(3) model, by means of the classical representation 
of the VBS stattP. 



(fik G [0, 2tt), we can define the the spin coherent state 
Ififc) as 



\n k ) = K4+^l) 2Sfc |0fe) 



(9) 



(|0fe) being the vacuum state at site fe), this states are 
complete but not orthogonal. Inserting the resolution of 
the identity 



2Sh 



1 



l 2S k + l 



4tt 



dn k \n k )(n k 



in M 



[A] 

{"}{/?} 



, and using the result 



-,Sk—liSk+l 



(10) 



n k ) = 



y/ (2S k )\u k h l v^ k+ \ the following form of the overlap 
matrix is obtained (dropping overall constant factors) 



M 



[A] 

MOT ~ 



i£A {i,j)EE A 
kGdA 

here fifc = (sin^fc cos tp k , smO k sin (p k , cos 9 k ) is the unit 
vector over the two dimensional sphere S 2 and u* is the 
complex conjugate of u. From[TT]we see that the overlap 
matrix M^l is hermitian, so the partial density matrix 
PA = Af(M^)*M^ is also hermitian. Using now that 
(wfe) 1_Qfc {v k ) ak = 0fe fc (abusing notation, a k goes from 
being a power, to become a (supra)index, ^ = u k , <f>\ — 
v k ) and the identity 



2^{4>lf = 5 af3 +n k -a a p, 



(12) 



where S a p is the Kronecker delta symbol, and a — 
(ai, CT2) 03) is a vector of Pauli matrices (no distinction is 
made between upper or lower greek indices); the expres- 
sion for the overlap matrix can be written as 



M 



[A] 

MOT 



/nfr n (i-^-^o as) 

ieA {i,j)eE A 

fcedA 

combining this result with Q , the density matrix of the 
subsystem A becomes 



IV. QUANTUM TO CLASSICAL MAPPING 

Introducing the spinor coordinates 4> k — (u k ,v k ) = 
(e^*/ 2 cos ^, e"^*/ 2 sin ^) at site k, with 6 E [0,2?r], 



fceA (i,j)eE A \jE B 



(14) 



<fc,i>e9 
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with Z the proper normalization factor to make TvpA = 
1. We can expand the matrix product inside (14) using 



the product identity for Pauli matrices Oi(jj = S^l + 
itijk&k (repeated index implies sum) where e^fc is the 
totally antisymmetric Levi-Civita tensor. The result of 
the term inside the square bracket in ( 14 1 is then 



[(I + 4'5)(I + nr?)W (15) 
= (1 + fij. • Cli)S af 3 + + £li + i(Cl k x Qi)) ■ a a p, 

where a x b represent the cross product between vectors 
a and b. 



V. 



GRAPH EXPANSION OF THE DENSITY 
MATRIX 



From (15), two types of expressions can be assigned 
to each edge on d. We draw an straight line between 
k and I whenever in that bond we have the expression 
(1 + Ofc - rii)6 a p, while we put a wiggly line for (fl, + 

fli + i(£lk x ^i)) ' &ai3- Expanding the product over the 
boundary in (|14|), we obtain a sum where each term has 



either a wiggly or straight line corresponding to (k, I) £ 
d. All the other bonds who don't belong to d have an 
straight line associated with them. 

In general, for a planar graph £, the expression 



j O{N) 



nf^ n (i+a*v%) (i6) 

kec N (i,j)&E c 



where VL is an N dimensional unit vector; corresponds to 
the partition function over C of the 0(N) mode. 27 which 
is analogous to a model of overlapping loops. To see this, 
we use that 



S 7 



N 



Sn 



odd 



■ n fc ) = o, (17) 



where Sn is the area of the S N ~ X sphere and the second 
property follows form the invariance of the integration 
measure under change £l k — ¥ — Cl,. As the only terms 
that contribute to Zq(n) are the ones with a product of 
even (f2j • O^) terms at each site of the graph, the whole 
partition function can be written as^ 



$>(C>C)- r(C) 



(18) 



with C a particular configuration of loops of total length 
T(C) that can be embedded in the graph C, and w(l,C) 
being the corresponding weight associated with a loop 
£ and with the particular configuration of loops C. For 
example, for the hexagonal lattice (coordination number 



Zi = 3) each site has associated just two bonds, and each 
integration of a site gives a factor of , except the last in- 
tegration which closes the loop. The partition function is 

then Z 0(N) = J2 C (f) 1 ^ Nn(C) > with n ( c ) the number 
of loops in the configuration C. The computation of spin 
correlations (f2 OT • O/.) corresponds then to the computa- 
tion of Zq(n)i with configurations that allow loops and 
open paths that begin at site m and end at site k. From 
(14) expanding the product over the partition's bound- 



ary we get a sum over different configurations of loops 
and open strands in the O(N) model, over the graph C 
with defects (wiggly lines). In the present case, x = — 1 
and N — 3 for the classical partition function of the VBS 
ground state. 

So far we have developed our ideas for general pla- 
nar graphs with no loops and no more than one bond 
shared between neighbors (Mij = 1), but from now on we 
will focus the discussion on translation invariant lattices 
with the previous restrictions. The discussion will remain 
general for lattices subject to the mentioned restrictions, 
that can be embedded on a torus. Using translation sym- 
metry, we can expand the product over the boundary in 
( 14 ) in different contributions of translational invariant 



Hamiltonians along the boundary, with increasing num- 
ber of non-trivial operators (Pauli matrices) acting on the 
local Hilbert space associated with a bond. The first term 
of the expansion correspond to the identity in the 2' a '- 
dimensional Hilbert space of the boundary. The second 
term, which is proportional to a constant external mag- 
netic field acting on the boundary chain, vanish. This 
follows from the observation that in this term, we have 
just one wiggly bond placed in the boundary - let's say 
at bond k with endpoints k and k - and the rest are 
just straight lines, which after integration will generate 
all the configurations of loops, and open lines that start 
at k, travel through the lattice and end at site k (for 
this type of bonds we will use dashed lines, to indicate 
the corresponding connection on the lattice). So we will 
have a term which is proportional to the spin correlation 
between k and k, and an integral of the form (see fig[3]a) 



dilk dClj. 

47T 47T 



(Sk ■ n k ) m (n k + Ci- k + tin, x %)) • s akh 



with m odd, which vanish trivially. The next terms in the 
expansion have two Pauli matrices acting on the differ- 
ent bonds, this terms are proportional to the only SU(2) 
invariants that can be constructed with two vectors (of 
Pauli matrices), namely Oi ■ <fj (see figwlb). Depend- 
ing on the separation between the wiggly bonds along 
the boundary, we have different contributions for which 
the numerical factor should decay exponentially with this 
distance, given that the VBS model is expected to have 
a mass jap (fact that is proven for linear and hexagonal 
lattices^) , result which is in agreement with the O(N) 
model being noncritical for N > 2 at x — — 1^. 
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FIG. 3. First terms in the graph expansion of the RDM pA- 
For a generic lattice, the number of bonds arriving to a bound- 
ary vertex (big circles) can be any even integer, here we show 
for simplicity the case corresponding to an hexagonal lattice 
where the number of bonds arriving to a vertex is exactly 2. 
Dashed lines show the remaining bonds after taking the trace 
over the whole lattice, except for the boundary vertices joined 
by wiggly lines. 



With the previous results, we can write the following 
expansion for the density matrix pa 



PA 



2 d 



-^A l]k a l -{S :j xS k ) + . . . (19) 



where I is the 2 d x 2 d identity operator and the coeffi- 
cients A r and are related to the correlation functions 
of the O(N) on the lattice C, with some sites and bonds 
erased along the boundary. Specifically for the first coef- 
ficient A r we have 




FIG. 4. (color online) a. Loop contribution to A\ in the 
hexagonal lattice. Big circles represent boundary sites, along 
the partition (dashed line), b. Dashed line represent a par- 
tition in the square lattice. The contributions to A\ consist 
now of configurations of overlapping loops. 



systems, we expect that this correlation decays expo- 
nentially with the separation of the spins, then we have 
A r ~ cxp(— r/^i). Numerical studies for two-leg VBS 
ladders have been performed^ being in agreement with 
this general result. For r = 1, the second term in (20) 



vanishes in the thermodynamic limit when minimum dis- 
tance paths joining the sites are cycles who travel the 
lattice. Also taking the limit of infinite size of the A and 
B subsystems, the interaction between the two bound- 
ary chains along different cuts of the partition vanishes. 
Then the total density matrix is the tensor product of 
matrices with the expansion (19), for each cut. 



It is clear that the first nontrivial term in the expan- 
sion ( [19] ) is the XXX Heisenberg Hamiltonian. We can 
also determine whether this interaction is ferro or anti- 
ferromagnetic in the simplest hexagonal lattice, from the 
loop expansion. The structure of the lattice determines 
the sign of the interaction through the number of bonds 
that define the allowed paths between site k and site fc + 1 
(each bond has an associated x — —1). An overall mi- 
nus sign comes from the contraction of two wiggly lines. 
Then it is easy to show that for the hexagonal lattice with 
a partition like the one in Fig. [4] a, all the paths connect- 
ing the boundary sites have even number of bonds, then 
the sign of the boundary XXX Hamiltonian is —1, so the 
boundary chain interaction is ferromagnetic. Numerical 
resultaSESlin finite size square lattices for a partition like 
Fig|4jb, indicate that in the square grid the interaction 
is anti-ferromagnetic. 



A r ~ ((r2 fc -r2 fc+r )(r!fe-r!fc +r )) £fc -((r! fe -r!fc +r )(%-r! fc+r )) £fc . 

(20) 

Here the correlation function is computed over the lat- 
tice C k which is the same lattice C but with the bonds 
(k, k) and (k + r, k + r) erased. This relation is exact 
for hexagonal lattices, while for other lattices with co- 
ordination number greater than 3, all the other possible 
contractions between even number of legs at the bound- 
ary sites have to be included. As usual with gapped 



Given the structure (19) of pA, we expect that the 
partial density matrix to contain all the commuting op- 
erators associated with the Heisenberg chain, which are 
also SU(2) invariant local operators, along with all their 
powers, in the form 



p A ~ exp 



d 

k=l 



(C)Q 



(fe) 



(21) 
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where the operators Q^> are obtained iteratively by the 
relation [B,Q (fe) ] = Q^ +1 \ (k > 2) with B the boost 
operator defined a d 30 * 31 * 



-Y 

2 ^ 



n\a r , 



>n+l 



1), (22) 



and = — \ Xm(^« ' &n+i — l)j the usual XXX Hamil 



tonian. While the relation (21 1 is just approximate valid 
at finite size, in the double scaling limit (d — > oo) we ex- 
pect it to be exact, being the coefficients of the expansion 
( 19 1 solutions of linear equations in terms of the IVm co- 
efficients. In this double scaling limit, pa should be not 
only SU(2) invariant, but also invariant under the full 
action of the Yangian V(S77(2))P. 



VI. CONCLUSIONS 

Given the structure of the VBS ground state, it is pos- 
sible to define on any planar graph, without loops, a VBS 
state, where the local spin at site i is given by Zi/2, with 
Zi the coordination number at site i. Using the Schwinger 
boson representation of the VBS ground state and the 
classical variable representation of this state, an expres- 
sion for the partial density matrix pa, which describe 
the physical subsystem A, obtained by partitioning the 
whole unique ground state, can be written. This expres- 
sion for pa decompose into a classical loop expansion 
of the 0(3) model in the gapped phase. Analyzing the 
different loop contributions, and assuming translation in- 
variance, we have shown that the partial density matrix 
that describes a subsystem of the VBS ground state can 
be expressed as a sum over different rotation-invariant 
quantum operators, where the Heisenberg interaction be- 



tween nearest neighbors gives the largest nontrivial con- 
tribution to the expansion. This quantum operators act 
on a spin 1 /2 chain in the boundary of the partition. The 
translational invariance assures us that the different con- 
tributions along the boundary are equally weighted, so 
the boundary operator is given by the XXX Heisenberg 
Hamiltonian. Here we discuss the case of translation in- 
variant lattices which can be embedded in a torus, but for 
other lattices with different topologies we expect similar 
results. 

For non translational invariant lattices, the first non- 
trivial local interaction term is expected to be also of the 
type ovdi+i but the Hamiltonian along the boundary will 
have different numerical prefactors for each local Heisen- 
berg interactions, generating a non invariant Heisenberg 
Hamiltonian in the boundary. 

We conjectured the possibility for the boundary chain 
to have an expansion in term of commuting charges of the 
XXX Hamiltonian and Yangian symmetry in the limit of 
infinite sites. 

The analysis shown in this paper should be useful for 
studying other dimensions d > 2 or other two dimen- 
sional lattices with more than one bond between a pair 
of sites. In that case, the local dimension of the spin 
operators in the boundary Hamiltonian should increase, 
having then boundary chains with higher representations 
of SU(2) per site, but still with SU{2) invariant local in- 
teractions. 

Another open problem is to see if this bulk-boundary 
correspondence extends for other points in the Hamilto- 
nian parameter space, close to the AKLT point. 
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